THE  PLATE  OVERLAP  TECHNIQUE 


/..  G.  TAFF 
Group  ')4 


ABSTRACT 


This  report  presents  the  development  of  the  plate  overlap  technique 
originally  suggested  by  H.  Eichhorn.  Included  are  discussions  of  standard 
coordinates,  plate  modeling,  the  plate  overlap  technique  itself,  and 
some  computational  details. 
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INTRODUCTION 


I . 

Over  the  past  three  years  we  have  been  systematically  bringing  together 
the  tried  and  true  techniques  of  classical  photographic  astrometry  and  the 
technology  of  today.  Our  aim  has  been  the  ultimate  amalgam  of  the  two.  The 
result  would  be  a completely  automatic,  real-time,  accurate  (-  1")  reduction 
of  the  position  of  any  celestial  object.  The  fulfillment  of  this  endeavor 
involves  the  work  of  many  people  and  the  supplementary  genesis  of  analytical 
tools,  complicated  software,  and  new  hardware.  Since  this  aim  was  clearly 
recognized  long  ago,  wherever  possible  in  earlier  documents  or  software 
development  flexibility  was  built  in.  Only  the  most  germane  parts  of  this 
material  will  be  repeated  here.  For  discussions  of  multi-color  photometry 
and  multi-color  systems  the  reader  should  see  1,  2,  and  3.  The  procedures 
one  should  use  to  update  stellar  positions  are  discussed  in  4 with  applica- 
tions to  the  SAOC  in  5.  Non-stellar  objects  are  discussed  in  6.  All  of  the 
pertinent  aspects  of  astrometry,  artificial  satellite  reductions,  elementary 
plate  modeling,  and  photographic  astrometry  are  in  7 and  8. 

The  main  purpose  of  this  report  is  to  present  an  even  more  powerful 
technique  of  data  reduction  which  is  known  as  the  plate  overlap  technique.  It 
was  originally  proposed  by  Eichhorn  (9).  His  idea  was  to  eliminate,  as  much 
as  possible,  systematic  errors  in  the  reduction  of  the  positions,  proper 
motions,  and  annual  parallaxes  of  stars.  While  the  conceptual  idea  is  simple, 
its  analytical  formulation  and  computational  implementation  are  not.  We 
discuss  this  in  §IV.  First  we  establish  our  notation  (§11)  and  review  the 
properties  of  a gnomonic  projection  (§111). 
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IT.  NOTATION 


We  use  the  standard  notation  of  the  subject  (unlike  in  8). 


a1?.,  »,  a™.  = plate  constants  in  the  xi  standard  coordinate 
i i kx  tj 

h>ijkf’  b"j  = plate  constants  in  the  eta  standard  coordinate 
c = the  color  of  the  m'th  star 


m 

hm  = the  number  of  epochs  for  which  a catalogue  position  of  the 
m'th  star  is  available 

llm  = the  total  number  of  catalogue  positions  for  the  m'th  star 

I,J  = upper  limits  to  indices  i and  j on  the  plate  constants 

m = the  magnitude  of  the  m'th  star 
m 

M = the  total  number  of  different  stars  on  the  n'th  plate 
n 

M = the  total  number  of  different  stars 

N = the  total  number  of  plates  on  which  star  number  m appears 

N = the  total  number  of  plates 

P see  F.q . (10) 

mn 

R see  Eq.  (11) 

m 

Rplate’  Rref  = t*le  residuals  from  the  plate  modeling  and  catalogue 
mismatch  of  the  reference  stars 


R = R 

t,  m = 
h 

t = 
n 

x 

mn 


ot,  6 
A,  A 


1 T I\  _ 

plate  ref 
the  epoch  numbered  h™ 

the  epoch  time  of  the  n'th  plate 

= the  measured  (not  necessarily  rectangular)  coordinates 

of  the  m'th  star  on  the  n'th  plate 

= the  equatorial  coordinates  of  a celestial  object 

= the  equatorial  coordinates  of  the  tangential  point 
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£,n  = the  standard  coordinates  related  to  a and  6 
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TIT.  THE  GNOMON TC  PROJECTION 

The  Imaging  system  (telescope  plus  camera)  projects  a portion  of 
the  celestial  sphere  onto  the  focal  plane  through  a point  on  the  optical 
axis  of  the  imaging  system.  By  analyzing  the  geometry  of  sucli  a projection 
we  can  go  from  coordinates  on  the  sky  (right  ascension  a and  declination 
A)  to  coordinates  on  the  focal  plane  (called  standard  coordinates  and 
universally  denoted  by  £,n).  When  the  unit  of  length  is  the  focal 


length  of  the  telescope  we  can  write 

cot6sin(a  - A)  = r, secA/ (n+  tanA) , (la) 

cot6cos(a  - A)  = (1  - r|tanA)/(n  + tanA),  (lb) 

tan(a  - A)  = £secA/(l  - ntanA) , (lc) 

£ ’ “ cotAsin(a  - A)/[sinA  + cot AcosAcos (a  -A)], 

= cosqtan(a  - A)sec(q  - A),  (2a) 

ri=Icos  A-  cot&inAcos(a  - A)j/[sinA  + cot  AcosAcos  (a  -A)], 

= tan (q  - A)  , (2b) 

cotq  = cot6cos(a  - A),  (2c) 


where  (A, A)  are  the  equatorial  coordinates  of  the  extension  of  the  optical  axis 
of  the  imaging  system  on  the  celestial  sphere.  This  is  also  known  as  the  point 
of  tangency.  The  angle  q is  a convenient  auxiliary  variable.  Also,  if 
|a  - A|  is  small,  one  should  use  Eq.  (Id)  in  place  of  Eqs.  (la,  lb); 

sinA  = (sinA  + ncosA)/(l  + ^ + r)^)^^.  (Id) 

The  orientation  of  the  ^,n  axes  is  such  that  the  positive  n axis  points 
north  and  the  positive  £ axis  points  east.  Equations  (1,  2)  are  completely 

•k 

See  SIVA  of  7 for  more  detail. 
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rigorous,  allow  one  to  speak  of  an  object's  position  In  either  coordinate 
system,  and,  just  as  the  equatorial  coordinates  are  measurable,  so  too  are 
tiie  standard  coordinates. 


5 


TV.  THE  PLATE  OVERLAP  TECHNIQUE 
A.  Motivation 

■k 

Suppose  we  have  several  plates  of  the  same  area  of  the  celestial 
sphere.  On  it  will  be  the  object  of  prime  interest  (the  program  object),  stars 
with  well  known  positions  and  proper  motions  (reference  stars),  many  other 
stars  (field  stars),  and  images  of  various  non-stellar  objects  (asteroids, 
planets,  galaxies,  etc.).  Let  us  suppose  that  each  plate  is  separately 


a model  i 

of  the 
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k=0 

1 

1=0 
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1,  2, 

3,.. 

n ..i  j k i 
a . . . „F,  n m c , 
lj  kf  mn  mn  m m 


, n j „k  t 

b . „£  n m c , 

ijkx  mn  mn  mm 


(3a) 


(3b) 


where  the  index  m = 1,  2,  3,...,  M numbers  the  stars,  the  index  n = 1,  2,  3, 

....  N numbers  the  plates  (taken  at  epochs  ft  }),  fx  ,y  1 are  the  measured 
r r n mn  mn 

coordinates  (not  necessarily  rectangular)  of  star  number  m on  plate  n, 

{ C , n ) are  the  computed  standard  coordinates  for  star  number  m at  epoch  t , 
mn  mn  n 

m is  the  magnitude  of  star  number  m,  c is  the  color  of  star  number  m,  and 
n m 

fa” jk£ ’^i jk£  ^ are  Plate  constants.  Note  that,  except  for  the  weak 
assumptions  of  a polynomial  form  for  the  model,  Eqs.  (3)  are  completely 
general . 

Once  we  obtain,  by  an  adjustment  procedure,  estimates  for  the  plate 
constants,  say  we  would  take  the  measured  values  for  the 


By  "plate"  we  mean  any  completely  formed  image  of  a portion  of  the 
celestial  sphere.  Whether  it  be  a photographic  plate  or  television 
camera  frame  is  Immaterial. 
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f f 

program  object,  (x  , } , invert  the  model  and  obtain  estimates  for  its 

I t 

standard  coordinates,  {£  FinaHy>  by  inverting  the  gnomonic  projection 

t I 

we  obtain  estimates  for  its  equatorial  coordinates,  {a  6 }.  If  only 

n n 

the  reference  stars  were  used  in  the  adjustment  for  the  plate  constants  (1)  we 
would  learn  nothing  more  about  these  stars,  (but  we  could  compute  the  positions 
of  all  of  the  field  stars  and  non-stellar  objects),  (2)  we  would  more  likely 
than  not  be  extrapolating  the  magnitude  and  color  terms  rather  than  inter- 
polating them,  ar.-d  (3)  if  the  time  scale  over  which  the  plates  were  taken  was 
sufficiently  small  to  render  the  effects  of  proper  motion,  foreshortening,  and 
parallax  factors  totally  insensible,  we  would  find  a total  of  N different 
values  for  the  position  of  the  program  object.  At  most,  one  of  these  would 
be  correct  and  the  remainder  would  be  systemat ically  in  error.  The  use  of 
averaging,  viz. 


N 

1 1 

a = 7 a /n, 

u n 
n=l 


£ = 7 6 / N , 

L n 


n=l 


(4) 


would  reduce  random  errors  by  a factor  of  l/\/u  but  not  the  systematic  ones. 

As  this  point  is  one  of  the  essential  reasons  for  inventing  the  plate 
overlap  technique,  a few  more  words  of  explanation  are  necessary:  Kven  if 
the  geometry  of  the  imaging  system  is  known  precisely  and  the  errors  of  the 
measuring  machine  modeled  exactly  [so  that  Eqs.  (3)  contain  both  the  necessary 
and  sufficient  number  of  parameters  and  an  incorrect  form  of  the  model  does 
not  systematically  bias  the  results],  Eqs.  (3)  do  not  hold  exactly.  The 
unavoidable,  random,  errors  of  measurement  imply  that  the  left-hand  sides  of 
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Eqs.  (1)  shou’d  be  written  as 


r 


t 


x - X 
mil  mn 


y - Y 
mn  mn 

where  X (Y  ) is  the  noise  associated  with  the  measurement  of  x (y 

mn  mn  mn  umn 

Since  a non-zero  mean  for  X (Y  ) can  be  handled  by  a redefinition  of 

mn  mn 

a11  (b'1  ) , we  assume  the  noise  to  have  a mean  of  zero  and  to  be 

oooo  oooo 


uncorrelated.  Hence,  our  adjustment  procedure  would  seek  to  minimize 


(5a) 


(5b) 

). 


R 

P 


1 ate 


N 

l 


n=l 


M 

v 

l 

m=l 


(X 

mn 


+ Y ) , 
mn 


(6) 


as  a function  of  fa"Jk£,  b^}. 

The  point  is,  because  of  the  existence  of  {X  , Y },  even  if  each  plate 

mn  mn  r 

images  the  same  area  of  the  celestial  sphere,  with  the  same  tangential  point, 
made  on  the  same  imaging  system,  measured  with  the  same  measuring  machine,  etc. 
it  will  not  be  true  that 


1 i j kf 


= A 


B 


= B. 


ijkT  "ijk*  " i.jkf 


(7) 


and  these  two  effects  will  not  cancel  when  computing  {£n,n  }.  Thus,  while  the 
program  object  clearly  has  one  and  only  one  position  on  the  celestial  sphere, 

we  can  not  deduce  it  and  each  set  of  computed  standard  coordinates  is 

» » 

systematically  off  from  (£  ,n  ). 

Now  consider  the  case  when  we  regard  not  only  the  plate  constants  to  be 
unknown  but  also  the  positions  and  proper  motions  of  al 1 of  the  celestial 
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objects  imaged  on  the  plate.  Moreover,  we  constrain  the  equations  of 
condition  by  insisting  that  the  reference  stars  have  their  catalogue 
positions  {f'  m,  n* , m}  at  the  epochs  ft.m),  h™  = 1,  2,  3,  Hm  for 

which  the  catalogues  exist.  Intuitively  we  feel  that  this  "bonding"  of 
the  plates  will  eliminate  all  of  the  systematic  effects.  Moreover,  the 
positions  and  proper  motions  of  the  reference  stars  are  automatically 
computed  and  the  color  and  magnitude  terms  are  not  extrapolated  since  all 
of  the  objects  contribute  to  the  solution.  Obviously  the  number  of 
unknowns  has  increased  dramatically,  the  analytical  and  computational  problem 
is  exponentially  more  complex  and  time  comsuming  (and  still  can't  be  done 
rigorously!),  and  our  ultimate  gains  may  be  small. 

When  the  plates  overlap  in  the  center-edge  pattern  (Fig.  la)  or  center- 
corner  pattern  (Fig.  lb)  the  power  of  the  technique  is  more  apparent.  Now 
fewer  reference  stars  are  needed  per  square  degree  to  obtain  the  same 
accuracy  since,  in  effect,  the  total  area  covered  by  the  overlapping  plates 
is  treated  as  one  large  plate. 

B.  Formulation 

In  order  to  proceed  with  a reasonable  economy  of  notation,  we  drop 

the  magnitude  and  color  terms  from  Eqs.  (3),  set  I =I'=I,J  =J'=J 

n n n n 

\fae[l,N),  and  ignore  the  possibility  that,  for  some  function  f,  a”.  = f({b”,  J-). 

The  first  assumption  is  only  to  assist  in  the  clarity  of  the  presentation.  The 

second  is  not  at  all  restrictive  if  we  regard  I = max(I  ,1'),  J = max(J  ,J') 

n n n n n n 

and  any  plate  constants  artificially  introduced  by  this  to  be  identically  zero 
and  not  subject  to  adjustment.  The  latter  restriction  is  substantiative 
(consider  the  four  constant  plate  model 
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x = aK  + hr]  + c. 


y = -bT  + art  + d 


(9a) 

(9b) 


which  allows  for  an  origin  shift,  a rotation,  and  a scale  change)  hut 

necessary  in  order  to  make  formal  analytical  progress.  We  shall  also 

assume  that  the  standard  coordinates  for  each  plate  are  referred  to  the 

same  tangential  point  and  let  h ' be  the  proper  motion  for  the  m’th  star 

m 

n 

in  the  f direction  and  |i  he  the  corresponding  quantity  in  the  n direction. 

As  before  we  will  have  a total  of  N different  plates  and  M different 

stars.  The  plates  are  numbered  by  n, the  stars  by  m,  and  the  n'tii  plate 

was  taken  at  epoch  t relative  to  an  arbitrary  epoch  of  t = 0.  On  the 

n't  1 1 plate  only  M different  stars  appear  (field  stars  plus  reference 

stars:  we  ignore  non-stellar  objects  in  the  sequel  as  they  add  nothing 

to  the  discussion).  The  total  number  of  different  plates  the  m'tli  star 

appears  on  is  denoted  by  . Finally,  to  allow  unrestricted  summations, 

we  define  I’  and  R by 
mn  m 


1 If  star  //m  appears  on  plate  II n 
0 otherwise 

1 if  star  II m is  a reference  star 
0 otherwise 


(10) 


(ID 


The  plate  model  is 
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f = f,  + M t = J y a..x  + " , 

mn  m m n /•  . ■ lj  mn  inn  mn 

1=0  ]=0 


(12a) 


n = r + j t = y ) b . . x y + H , 

mn  m m n „ tj  mn  mn  mn 

1=0  i=0 


(12b) 


Vm^fl.M],  nc[l,N]  3 P ^ 0. 

mn 


The  reversal  with  respect  to  Eqs.  (3)  and  the  precise  meaning  of  the 

noise,  {n  ,11  ) is  discussed  in  § IVC  below.  That  part  of  the  residuals 

mn  mn 


due  to  the  plate  modeling  is 


N M 

r , = y y p [(-  ) + (h  r]. 

plate  ' L mn  mn  ran 

n=l  m=l 


(13a) 


The  other  part  of  the  residuals  are  from  the  discrepancy  between 

c c 

the  catalogue  position  of  a reference  star  number  m,  (£  ,m,  n , m)  at  epoch 

mh  mh 

t = t^  m for  a total  of  lim  different  catalogue  positions,  e.g., 


K + l/'t  m = F,C  m + 7iC.m, 
m m h mh  mh 


(14a) 


n c c 

n + V t.  m = n t in  4-  H , m, 
m m h mh  mh 


(14b) 


hm  = 1,  2,  3,  ....  H1"  Vmr.fl.M]  3R  t 0. 


Hence,  that  part  of  the  residuals  due  to  the  catalogue  positions  is 
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(1'Jb) 


R = )’  H [c^m)2  + OlC.  m)5!, 

ref  ( ml  mil  mil 


R = R , + R r 

p 1 a t e re  f , 


(13c) 


If  R = f)  then  h"  llm  0 and  C . in,  ri*  m are 
m mn  mli 


x and  y are  undefined, 
mn  mn 


The  normal  equations  are  given  by 


unde f i ned . 


Similar  1 y , if  P 

mn 


3R/3a"  = 0 

i 1 

3R/hbn,  = 0 

i i 


n = 1 , 2,  3,  . . . , N, 


(15) 


i = 0,  L,  2,  ....  I,  j = 0,  1,  2,  ....  J; 


3R/3T  = 0 

m 

<)  R / c)r  I = 0 

m 

'jR/^j  * 0 
m 

3R/3(jn  = 0 
m 


m = 1,  2,  3,  .*•>  M. 


(lb) 


II 


,m 


With  the  understanding  that  ) f (hm)  = 0 if  R =0,  these  equations 

m 

■ m . 

it  ii  =1 

ran  ht*  written  as 


M 

y p„ 

m=l 

M 

y p 

i n 

m=l 


r e n i+k  S+t 
l_0  ai.iXn,n  ymn 


I .1 

y i 

i =0  j=0 


l 1.  y, 


n i+k  i+f 
i)Xmn  ymn 


M „ 

- y P £ x y = 0,  (17a) 

'•  mn  mn  mn  mn 
m=l 


M . t 

- } P n x y =0,  (18a) 

' mn  mn  mn  mn 
m=l 


n = 1,  2,  3,  ....  N,  k = 0,  1,  2,  ....  I, 

l = 0,  1,  2 1; 
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N 


T J 


y p y y an.x*  yJ  + r <n  + nm)  + 

mn.1-  . „ li  mn  mn  m m 


n=  1 t=0  j=0 

N 


r,  ¥ I!"1  \ H" 

1Jm  'mn1.,  + A V/  = ^ W’ 

V‘=1  hm=l  ' hm=i 


(17b) 


) I / / h..x  y + 0 (N  + H ) + 

‘ , mn.1'  . • i]  mn  mn  m m 

1=1  1=0  j=0 


/ N Mm  \ Hm 

li  1 1 } I*  t + y t m = y n\  m, 

\n=1  hm=l  ' hm=l 


(18b) 


N 

- I 

n=l 


I .1 


- P t / ) a . . x y + P t + y t 

, mn  n . ' ii  mn  mn  ml  , mn  n ' 1 

n=l  1=0  i=0  \n=l  , m . 

' h =1 


N 


H 


H 


+ U y P t + y tm|=y  £ , mt  m, 
m , mn  n hi  '■  mh  h 

n=l  , m , / . m 

h =1  ' h =1 


I J 


N 


H 


y p 1 y y b^.x*  + n [ y p t + y t 

u mn  n . _ . ' ii  mn  mn  ml  , mn  n u r 
i=0  j=0  \n=l  hm=1 

+ y 1 f / P t + y t m ] = y r|c  mt  m, 
m I '■  mn  n L h L mh  h 

\n=1  l,m=l  / hm=l 


(17c) 


(18c) 


in  = 1 , 2 , 3,  . . . , M . 

Th  i s completes  die  analytical  specification  of  the  problem. 

C.  On  Statistical  Rigor 

Suppose  two  qantities,  u and  v,  are  linearly  related,  say 

u = kv  . (19a) 
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Further  suppose  we  have  I1  sets  of  values  of  u and  v,  (u  , v } with  random 

P P 

errors  i in  u . Hence,  the  correct  equations  of  condition  are 
P P 


u =kv  + C , p = 1 , 2 , 3,  ....  P. 
P P P 


(20) 


We  determine  an  unbiased  estimator  for  k (K)  by  minimizing 


S = y r.2  = y (u  - kv  )2 
p=l  P p=l  P P 


(21a) 


with  respect  to  k.  The  result  is 


P P 

1C  = y u v / y v 

^ p=]  P P 


2 

q=l  q 


(21b) 


P P - 

= k + / v t /}  v . 

p p , q 

p= 1 q= 1 


(22a) 


When  performed  in  this  fashion  the  full  set  of  results  available  from 
least  squares  theory  is  available  to  us  and  complete  statistical  rigor  is 
satisfied.  In  particular,  we  know  that 


1 imK  = k . 

P 

KI 


(23a) 


Had  the  relationship  been  written  as 


v = <u 


(19b) 


with  u still  subject  to  random  error  (but  v not),  then  blindly  applying 
least  squares  procedures  to 
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(21b) 


/ (v  - KU  ) 
p=  I P [> 


2 


resu Its  in 


1/K  = )'  v u / }’  u 

|>  p p ' . q 

p=l  q= 1 1 


P 

= k 1 

P=1 


v C 
P P 


2 

v 


P 

+ 2k  l 

q=l 


v i 

q q 


p 

+ 7 

q=l 


2 

(, 

q 


(22b) 


.111(1 


lim  1/K  = k <v2>/[k2<v2>  + <l2>] 


4-  1/k. 


(23b) 


Hence,  the  estimator  for  K Is  biased. 

Thus,  F.qs.  (3,  5)  are  rigorous  and  F.qs.  (12)  are  not.  Moreover,  the 

exact  interpretation  of  {5  ,H  } Is  unclear.  Unfortunately  the  rigorous 

n mn 

equations  yield  a non-linear  estimation  problem  for  which  there  is  little 
theoretical  basis.  In  addition,  the  resulting  normal  equations  are 
horrendously  complicated.  We  therefore  have  the  choice  of  posing  (correctly) 
an  insoluble  problem  or  posing  (incorrectly)  a soluble  one.  We  have  chosen 
tilt;  one  which  is  computationally  easier  to  handle  (e.g.,  the  latter). 
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I). 


On  Solving  the  Normal  Equations 


It  would  be  cavalier  of  us  to  remark  that  the  normal  equations 
are  linear,  and,  therefore,  standard  techniques  are  applicable.  Since  their 
order  is  2 N ( I + 1)  ( I + 1 ) + AM  one  can  easily  have  thousands  of  unknowns. 
When  the  pattern  is  of  the  center-edge  or  center-corner  type  the  matrix  of 
the  equations  is  of  the  "banded-bordered"  type  and  special  techniques  are 
available.  If  one  is  willing  to  iterate  towards  a solution  then  the 

A 

problem  is  even  simpler.  In  Eqs.  (17b,  1.7c)  first  set  I'  and  N equal  to 

1 inn  m 

zero.  Solve  these  for  the  first  approximation  {£,  ,\i  ' } Vl<  /0 . We  now  solve 

mm  m 

Eq . (17a)  with  P / 0 iff  R ^ 0.  This  provides  the  first  approximation 
mn  m 

to  far'.}.  We  can  now  compute  {£  } for  the  field  stars  and  insert  these 

i j mn 

values  into  our  simplified  versions  of  Eqs.  (17b,  17c)  for  all  of  the  stars. 

When  this  system  is  solved  wo  have  obtained  our  second  approximation  to 

If  , u '}  for  all  of  the  stars.  We  then  return  to  solving  for  the  star 
m m 

constants  using  all  of  the  stars,  etc. 


« n r n 

The  assumption  that  a^.  ^ f(lb.'.t))  separates  f,  from  r)  and  we  only  consider 
the  f coordinate  here.  The  n coordinate  is  treated  in  an  analogous  fashion. 
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V.  ASTRONOMICAL  APPLICATIONS 

The  astrophysical  benefits  of  the  plate  overlap  technique  have  been 
noteworthy.  Gatewood  and  his  collaborators  at  Allegheny  Observatory 
have  rapidly  applied  the  technique  to  Allegheny's  extensive  plate  collection. 
Among  their  major  contributions  have  been  (1)  a complete  rediscussion  of 
the  hypothetical  planet  around  Barnard's  star,  (2)  a test  of  the  theory 
of  general  relativity  and  stellar  structure  theory  by  astrometrical ly 
inferring  the  gravitational  redshift  of  the  white  dwarf  van  Maanen  2, 
and  (3)  the  beginnings  of  a thorough  rediscussion  of  the  anomalous 
multiple  systems  used  to  define  the  zero-age  main  sequence  mass-luminosity 
relationship  (9  Pup,  r)  Cor  Bor,  and  10  UMa) . Perhaps  similar  gains  in 
our  knowledge  of  the  geopotential,  the  selenopotential , solar  radiation 
forces,  etc.  await  us. 
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